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ABSTRACT 

We develop a new method to predict the density associated with weak lensing maps of (un)relaxed 
clusters in a range of theories interpolating between GR and MOND (General Relativity and Modified 
Newtonian Dynamics). We apply it to fit the lensing map of the bullet merging cluster 1E0657-56, 
in order to constrain more robustly the nature and amount of collisionless matter in clusters beyond 
the usual assumption of spherical equilibrium (Pointecouteau & Silk 2005) and the validity of GR on 
cluster scales (Clowe et al. 2006). Strengthening the proposal of previous authors we show that the 
bullet cluster is dominated by a collisionless - most probably non-baryonic - component in GR as 
well as in MOND, a result consistent with the dynamics of many X-ray clusters. Our findings add 
to the number of known pathologies for a purely baryonic MOND, including its inability to fit the 
latest data from the Wilkinson Microwave Anisotropy Probe. A plausible resolution of all these issues 
and standard issues of Cold Dark Matter with galaxy rotation curves is the "marriage" of MOND 
with ordinary hot neutrinos of 2eV. This prediction is just within the GR-independent maximum of 
neutrino mass from current /3-decay experiments, and is falsifiable by the Karlsruhe Tritium Neutrino 
(KATRIN) experiment by 2009. Issues of consistency with strong lensing arcs and the large relative 
velocity of the two clusters comprising the bullet cluster are also addressed. 
Subject headings: gravitation - dark matter - galaxy clusters - gravitational lensing 



The bullet interacting cluster IE 0657-56 has recently 
been argued to have produced the first completely unam- 
biguous evidence that galaxy clusters are shrouded in a 
dominant component of collisionless dark matter (Clowe 
et al. 2006, hereafter C06). The explanation for this was 
that the peaks of the convergence map are offset, without 
a shadow of a doubt, from the main observable baryonic 
components, i.e., the gas marked by the bright X-ray 
emission. Instead, the lensing signal peaks at the galax- 
ies (the minor observed baryonic components, about 1/7 
of the total X-ray gas mass) which is exactly where we 
expect any collisionless dark matter (DM) halos to center 
on. It was also argued that modified gravity theories, in- 
cluding the relativistic counterpart of MOND (Milgrom 
1983, Bekenstein 2004), would have no way of reproduc- 
ing the lensing map since the gravity in these theories is 
thought to trace the light. 

There are, however, two caveats with this line of rea- 
soning. As first cautioned by Angus, Famaey & Zhao 
(2006, hereafter AFZ06), the features in the lensing con- 
vergence map in a non-linear gravity theory do not al- 
ways reflect features in the underlying matter surface 
density in highly non-spherical geometries. For exam- 
ple, in MOND, the convergence can indeed be non-zero 
where there is no projected matter (Zhao, Bacon, Tay- 
lor & Home 2006, Zhao & Qin 2006). The other caveat 
is whether ordinary collisonless neutrinos (which are de- 
tected particles known to have a small mass, although 
non-baryonic and non-photon-emitting as all leptons in 
the standard model arc) should be given equal status as 
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the "known" matter (e.g., the hadrons in gas and stars) 
or the stigma of "unknown" dark matter (e.g. primodial 
black holes, exotic cold WIMPs from SUperSYmmetry). 
Indeed, ordinary neutrinos represent at least 0.1-3 times 
the density of "known" gas and stars in a galaxy cluster 
given the current experimental mass range of A'^ > 3 fla- 
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As for MOND, the real question is actually to examine 
whether the bullet cluster poses any new challenges to 
MOND at galaxy cluster scales. It is indeed well-known 
that the dynamical mass from the X-ray temperature 
profiles of clusters in MOND exceeds their baryonic con- 
tent (Aguirre et al. 2001). As a fix, a component of 

1 — 2eV neutrinos has been invoked to explain cluster 
cores (Sanders 2003) and the Cosmic Microwave Back- 
ground (McGaugh 2004). 

Nevertheless, applying a new method developed to fit 
the detailed weak lensing map of the bullet cluster, we 
place robust constraints on the dark matter density in 
non-spherical non-equilibrium geometry. We test a range 
of gravity theories interpolating smoothly between MOND 
and GR. We call these models the /Lt-gravity Hot Dark 
Matter models (/xHDM, with appropriate interpolating 
function /i and hot neutrino content) in contrast to its 
contending ACDM model (with appropriate cosmologi- 
cal constant A and cold DM content). Note both cos- 
mologies can drive late-time acceleration of the universe 
(Diaz- Rivera et al. 2006, Zhao 2006), and form struc- 
tures and CMB anisotropics (Skordis et al. 2006, Dodel- 
son & Liguori 2006). 

GR is a limiting case of the multi-field theory of grav- 
ity, TeVeS (Bekenstein 2004) as Newtonian gravity is a 
Hmiting case of MOND (Bekenstein & Milgrom 1984, 
BM84 hereafter). In these theories, the total potential 
is due to the usual Newtonian potential of baryons plus 
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a baryon-tracking scalar field (sec e.g. Zhao & Famaey 
2006), which creates the DM or MOND effect. While 
there are subtle differences with MOND in non-spherical 
geometries, in the limiting case of scale-free flattened 
models, AFZ06 showed that using the Poisson-like equa- 
tion of BM84 for the total gravitational potential $ was 
a reasonable approximation to the multi-field approach. 

We thus investigate here a class of MOND-like non- 
linear laws of gravity, where the gravitational potential 
$ satisfies the Poisson-like equation of BM84. The aver- 
age mass density p(< r) or the total mass (e.g., baryons 
and neutrinos) of the system enclosed inside any radius 
r centered on any position will thus be estimated from 
the divergence theorem with an effective gravitational 
parameter Ggg (not a constant): 



Mbary (r ) -FM^(r) = J 
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where the surface area element dA = sm{9)d9dip and 
the interpolating function (cf. AFZ06) q^^(^^) = fJ-ix) = 
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case a ^ oo (i.e. Gcff = G) corresponds to General 
Relativity, whilst a = is the toy-model (Eq.64) of 
Bekenstein (2004), and a = 1 is the simple /U- function 
which has a better fit to the terminal velocity curve of 
the Milky Way (Famaey & Binney 2005). We also ex- 
amine the standard interpolating function /i(x) = 

for comparison with other works. 

Fitting the convergence map of a multi-centred X-ray 
cluster: In GR, the convergence map allows us to imme- 
diately derive the underlying projected density of matter. 
However, as shown in AFZ06, the situation is different 
in MOND, where what you see (in terms of convergence) 
is not always what you get (in terms of density). For 
that reason, we use a potential-density approach here- 
after: we fit the convergence map using a parametric set 
of potentials, and then use the best-fit potential to derive 
the corresponding surface density for various choices of 
the gravity's interpolating function fi. 

The bullet cluster is 4-centred, the centres being the 
positions of the main cluster's CoUisionless Matter (re- 
ferred to as CMl, including its member galaxies) and 
X-ray gas components (XRl), and the sub-cluster's CoUi- 
sionless Matter (CM2) and X-ray gas components (XR2). 
We chose to model those 4 mass components as 4 spheri- 
cal potentials: note however that, in non-linear gravities, 
the 4 mass densities corresponding to those spherical po- 
tentials will not linearly add up, especially when fj, is 
rapidly varying with position inside the system. 

We thus write the lens-potential as a superposition of 
four potentials pinpointed at four centres ri: 



$(f) 
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Each potential is fully described by two parameters, the 
asymptotic circular velocity Vi and the scale length pi. In 
GR these potentials correspond to cored isothermal den- 
sity profiles. These potentials are similar to the cusped 
potentials of AFZ06, which were found to have less good 
fit to the convergence map. 



Using Fig. lb of COG we set up a coordinate system for 
the bullet cluster. The centers XRl, XR2 and CM2 he, 
to a first approximation, along the RA direction, which 
we chose as our a;-axis. Our ^;-axis is along the line of 
sight. As suggested by Markevitch et al. (2004) and C06, 
we chose the four centres of the potential to be exactly 
in the x-y plane with their (AT, Y) coordinates chosen at 
the four observed peaks. As a consequence, the potential 
of Eq.(2) has 8 parameters. 

The parametric convergence map in the x — y plane is 
simply computed by linear superposition of the individ- 
ual contributions to the convergence from the four spheri- 
cal potentials (see AFZ06) , the convergence of each solely 
depending on its parameters vt and pi, and the rescaled 
radius Si=pr^^/{x-Xi)^-\-{y-Yi)^: \ 
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We then reproduced the observed convergence map 
(Fig. lb of C06) by least-squares fitting the asymptotic 

velocities Vi and concentration parameters pi of each of 
the spherical potentials, using n = 233 points from the 
K = 0.16, 0.23, 0.3, 0.37 contours (always with a constant 
number of points per contour length). We also tried mov- 
ing the 2-coordinate centres well within the errors of the 
brightest cluster galaxy. 

The goodness of fit of the model Kmodei.i to ob- 
servations Kobs,i is given by a characteristic variance 
if . Due care was taken to max- 
imize the resemblance to the X-ray gas features with the 
centres as marked in Fig. 1. Another constraint was 
trying to ensure a reasonable mass of X-ray gas to con- 
form with the estimates of C06 and Bradac ct al. (2006, 
hereafter B06; see discussion section). The best-fit pa- 
rameters (listed in Table 1) yield the convergence map 
shown in Fig.l upper panel, with a variance ^0.1 x0.07^, 
which is acceptable since 0.07 is both the typical obser- 
vational error at individual points and the convergence 
spacing between two neighbouring contours. 

Masses of gas and collisionless matter in various postu- 
lated gravities: Applying Eq. (1) to our potential model 
we can predict the matter volume density in the clus- 
ters, e.g., the values given in Table 2 up to a trivial 
degeneracy. Integrating over the line of sight, we note 
that the projected density contours are slightly different 
from that of convergence contours in non-linear gravities 
(cf. dashed blue contours of Fig.l). While confirming 
AFZ06, this non-linear effect appears much milder than 
expected earlier. 

In order to match the observed X-ray gas mass, which 
is a minor contributor to the Icnsing map, we use the 
asymmetry in the calculated surface density to subtract 
off all the collisionless matter centred on the galaxies 
(CMl and CM2). The key here is to notice the symme- 
try of galaxies around the dashed line joining the centres 
of the two galaxy clusters (cf. Fig. 1 upper panel) . If we 
fold the map over the axis of symmetry subtracting the 

^ Erratum of AFZ06: correcting a typo, their Eq.25 should read 



(\R-R,\^ -2pf)D,l!B{\R-Ri\) 



'^^^^ 2|i?-i?i|{|fl-Hi|2-p?) ^ c2(|fl-H,|2_p2) 

bending angle 6 is given by Eq.23 of AFZ06 
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TABLE 1 

Best fit parameters of the convergence map. Squared velocities in (1000 kms~^)^ and scale radii, p in kpc. 
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Fig. 1. — Our fitted convergence map (solid black lines) overplotted on the convergence map of C06 (dotted red lines) with x and y 
axes in kpc. The contours are from the outside 0.16,0.23,0.3 and 0.37. The centres of the four potentials we used are the red stars which 
are labelled. Also overplotted (blue dashed line) are two contours of surface density [4.8 & 7.2]x10^A^q pc"'^ for the MOND standard /i 
function; note slight distortions compared to the contours of k. The green shaded region is where matter density is above 1.8 X 10"'^ Mq pc~^ 
and correspond to the clustering of 2eV neutrinos. Inset: The surface density of the gas in the bullet cluster predicted by our coUisionless 
matter subtraction method for the standard /i-function. The contour levels are [30, 50, 80, 100, 200, 3OO]M0pc-2. The origin in RA and 
dec is [06''58™24.38^-55°56'.32] 



TABLE 2 

Compares the results of COO and B06 with estimates of projected mass in [10-'-'' Mq] for gas around the X-ray centres and 

TOTAL mass around THE LENSING CENTRES IN THREE DIFFERENT GRAVITIES (OR, STANDARD AND SIMPLE fj,); A BIAS FACTOR OF 1.15 — 1.73 
SHOULD BOOST OUR asymmetry-based GAS MASS ESTIMATES HERE. BOO USED AN ELLIPSE "WITH SEMI-AXES OF 250 KPC AND 150 KPC AROUND 

XRl, WHICH WE COMPARE WITH A CIRCLE OF 180 KPC (THE SAME AREA). THE LAST COLUMN PREDICTS THE AVERAGE MATTER DENSITIES 
IN [10"'^ M0pc~^] WITHIN IOOKPC of THE LENSING PEAKS; THE DEGENERACY OF LENSING TO A STRETCHING OF THE POTENTIAL ALONG THE 
LINE OF SIGHT 4>{x,y,z) (<f>{x,y, zd^) MEANS THAT THE CENTRAL DENSITY PREDICTIONS WOULD BE EASILY lowered IF WE ADOPTED C < 1- 
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r<100/180kpc 


r<100/80kpc 
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r < 100 
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GR 


1.05/1.97 


0.33/0.27 
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17.2 


2.63 


2.59 


standard /i 


0.97/1.79 


0.29/0.24 


9.0 


6.78 


2.26 


2.34 


simple /I 


0.74/1.33 


0.21/0.18 


7.13 


6.42 


1.66 


1.76 


C06/B06 


0.66/2.0 


0.58/0.42 


20.0/28.0 


21.0/23.0 







lower part from the upper part we are left with the ma- 
jority of the gas. Then we performed a straight forward 
numerical integration over the areas defined above. The 
lower panel of Fig.l demonstrates that this technique 
works well in separating the surface density of gas from 
the collisionless matter. The values for the gas mass for 
our three gravities are given in Table 2. Note this tech- 
nique works less well for the sub cluster (XR2) as it lies 
quite close to the axis of symmetry and thus much gas is 



cancelled out by other gas. For GR only we can directly 
compare the gas corresponding to the potential and that 
calculated by our subtraction method. For the main clus- 
ter, we find that integration of the surface density gives 
2.32 X IO^^Mq within the ISOkpc aperture which is 15% 
more gas than estimated by asymmetry. For the sub 
cluster we find 5.7 x lO^^M© from integration within the 
lOOkpc aperture of the gas center, 73% more than from 
asymmetry. As such, in MOND we can expect the gas 
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masses to increase by similar amounts and this helps to 
explain the low gas masses found, especially in the sub 
cluster (XR2). The reason our K-map is skewed towards 
the gas peaks is a feature of the cored isothermal po- 
tentials. Table 2 shows we pack too much gas into the 
central lOOkpc of the main cluster compared to that ob- 
served only for it to balance by ISOkpc. Using a potential 
that correctly matches the gas density would not skew 
the map. 

Table 2 also compares the BOG and COG projected mass 
within a 250 kpc circular aperture centred on both galaxy 
clusters with our total mass within these apertures for 
three gravities (GR, simple fi and standard /i). Clearly, 
these amounts of mass exceed the observed baryons in 
gas and galaxies over the same apertures, by a factor of 
3 even in MOND. While very dense clumps of cold gas or 
MACHOs are still allowed by the missing baryon budget 
to reside in galaxy clusters without much collisions, we 
will focus on the possibilities of fermionic particles being 
the unseen collisionless matter in the lensing peaks. 

Following Tremaine & Gunn (1979), we use the dens- 
est regions of the collisionless matter to set limits on the 
mass of ordinary/sterile neutrinos. A cluster core made 
of neutrinos of mass would have a maximum density 

(Sanders 2003) satisfying ^^^^^Smopc--^ (liv)"'^' = 

v^ov) "^tiere we adopted a mean temperature of 9 keV 
in the two clusters. Comparison with the regions of the 
highest volume density of matter shown in Fig.l upper 
panel suggests that the relatively diffuse phase space den- 
sity in the bullet cluster is still consistent with ordinary 
2 eV neutrinos making up the unseen collisionless com- 
ponent. Note the lensing- predicted 3D density is gener- 
ally non-unique due to the degeneracy of a line-of-sight 
stretching of the potential. While a better fit to the gas 
mass and the lensing map could be produced by using 
several ellipsoidal potential components as opposed to 
the rigid four spheres with fixed centers here, the present 
model suffices as a demonstration. 

Discussion The lensing-predicted 2 eV mass for neu- 
trinos in the non-equilibrium bullet cluster greatly 
strengthens the finding of Pointecouteau & Silk (2005) 
that (rrij/ > 1.6 eV using standard /i) based on spherical 
gas equilibrium of other clusters. 

Moreover, the 2eV neutrinos are falsifiable in the 
near future. At present it is compatible with model- 
independent experimental limits on electron neutrino 
mass m^^e < 2.2 eV from the Mainz/Troitsk experi- 
ments of counting the highest energy /3-decay electrons 

of — » '^He^ + e~ + i/e + 18.57 keV (the more massive 
the neutrinos, the lower the cutoff energy of electrons). 
The KATRIN experiment (under construction) will be 
able to falsify 2eV electron neutrinos at 95% confidence 
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within months of taking data in 2009. Our prediction of 
a 4th (hot sterile) neutrino is fasifiable by the Booster 
Neutrino Experiment. 

We note finally a couple of discontinuities with the 
work of COG and BOG. Our adopted lensing map of COG 
(with a peak k = 0.37) implies a surface density that 
is too weak to form the observed large scale arcs in any 
gravity for sources at any redshift. Secondly, our MON- 
Dian models greatly reduce the amount of collisionless 
matter needed to fit the map of BOG in GR; a reduction 
by 3-4 times at 250 kpc. These masses are all integrated 
over the line of sight and as such give poor estimates 
of the mass in the system for comparison with the grav- 
ity independent gas masses. Consistency between strong 
and weak lensing data remains to be understood together 
with issues of smoothing, normalization and zero point 
of the K maps. 

Nonetheless, the data still convincingly require a dom- 
inant component of collisionless and most probably non- 
baryonic matter at cluster scales. A traditional miscon- 
ception is that the existence of a large quantity of non- 
baryonic matter would make a modified gravity theory 
such as MOND contrived or redundant. This is not the 
case with ordinary hot neutrinos, which are too diffuse to 
either perturb the good MONDian fits to galaxy rotation 
curves nor explain these curves in GR. 

As a tie-breaker between ACDM and /iHDM we note 
that high-speed encounters are rare in CDM simulations 
(Hayashi & White 200G). The potential well of the 
main cluster would be too shallow to accelerate the sub- 
cluster without stronger-than-Newtonian gravity (Far- 
rar & Rosen 200G). Our MONDian isothermal poten- 
tial well (cf. Eq. 1 and Table 1) would accelerate 
the two clusters to a maximum relative speed Wmax = 
y^vf ln(l + (r/pi)2) + v'i ln(l + (r/pz)^) - 4800km/s ii 
^ 33v — ~ 1-8 Gyrs ago the two clusters turned 
around from the Hubble flow at r ~ 2rvh- ^ 2 x 1500 
kpc and free fall towards each other along the east-west 
direction. 
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